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Abstract 

We study theoretically the optical response of a nanohybrid comprising a symmetric quantum 
dimer emitter coupled to a metal nanoparticle (MNP). The interactions between the exitonic 
transitions in the dimer and the plasmons in the MNP lead to interesting effects in the composite’s 
input-output characteristics for the light intensity and the absorption spectrum, which we study in 
the linear and nonlinear regimes. We find that the exciton-plasmon hybridization leads to optical 
bistability and hysteresis for the one-exciton transition and enhancement of excitation for the two- 
exciton transition. The latter leads to a significant decrease of the field strength needed to saturate 
the system. In the linear regime, the absortion spectrum has a dispersive (Fano-like) line shape. 
The spectral position and shape of this spectrum depend on the detuning of the dimer’s one-exciton 
resonance relative to the plasmon resonance and the ratio of the exciton-plasmon coupling constant 
and the exciton dephasing rate. Upon increasing the applied field intensity to the nonlinear regime, 
the Fano-like singularities in the absorption spectra are smeared and they disappear due to the 
saturation of the dimer, which leads to the MNP dominating the spectrum. The above effects, for 
which we provide physical explanations, allow one to tailor the Fano-like shape of the absorption 
spectrum, by changing either the detuning or the input power. 

PACS numbers: 78.67.-n 73.20.Mf 85.35.-p 
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I. INTRODUCTION 


Over the past decade, heterogeneous nanosystems, consisting of quantum emitters (quan¬ 
tum dots, molecules, molecular aggregates, etc.) in close proximity to metal surfaces^*^ or 
metal nanoparticles (MNPs), have received a great deal of attention. The intriguing fea¬ 
tures of these systems arise from the hybridization of different types of optical excitations, 
in particular, excitons (in the quantum emitter) and plasmons (in the metal). The exciton- 
plasmon coupling can drastically modify the optical response of hybrids, leading to inter¬ 
esting physical phenomena, such as optical bistability,-^— exciton-induced transparency!^ 
enhancement of Rabi oscillations,— suppression of quantum coherence via infrared driven 
coherent exciton-plasmon coupling,—!!^ florescence,—*!^ Forster energy transfer,— photo¬ 
luminescence quenching,— Fano-like absorption,-*2ir— and other exciting effects.— These 
phenomena may have a strong impact on the development of active nanophotonic devices 
and metamaterials (e.g., optical switches, single photon sources, biosensors, etc.). 

In this paper, we theoretically investigate the optical response of a nanocomposite con¬ 
sisting of a symmetric dimer of two two-level molecules coupled to an MNP. This work is 
an extension of works performed earlier by us^*!^ and other authors,-*!^ on the optical re¬ 
sponse of single two-level molecules coupled to an MNP. These earlier studies revealed the 
interesting possibility of bistable optical response for such nanocomposites, arising from the 
self-action of the molecular excitation on itself through the reflection of the electric field 
generated by it on the MNP. This is a direct consequence of exciton-plasmon coupling. 

The dimer is a first step to considering molecular aggregates coupled to MNPs, a topic 
that recently has attracted experimental interestIn contrast to a single two-level 
molecule, a dimer represents an optical ladder system consisting of a ground state, two 
one-exciton states, and one two-exciton state. Direct (one-photon) transitions between the 
ground state and the two-exciton state are not allowed. By their richer excited-state struc¬ 
ture, dimers coupled to MNPs may be expected to reveal new effects in their optical response, 
as compared to single molecules coupled to MNPs. Examples inspired by effects observed 
for ladder type (three-level) quantum dots are two-photon Rabi oscillations^*^^ and suppres¬ 
sion of quantum decoherence.— i!^ Another effect that in principle could occur is multi-stable 
behavior.— 

We focus on two types of aspects of the optical response of the hybrid. The first one is the 
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intensity dependent response at a few selected driving frequencies. Quantities of interest are 
the occupation probabilities of various exciton levels in the system and the expectation value 
of the dipole of the dimer. We will show that the one-exciton manifold gives rise to bistable 
response and hysteresis in these quantities, while the two-exciton manifold does not. No 
multistability is found. We also show that two-photon absorption to the two-exciton state 
is strongly enhanced by the exciton-plasmon coupling, an effect that also has been observed 
experimentally^ and which is of interest to two-photon microscopy. This enhancement also 
leads to a low saturation intensity of the dimer as compared to the case where it is not 
coupled to a MNP. 

The second aspect we investigate is the absorption spectrum of the hybrid system. Most 
of the results obtained in this part do not depend on the presence of a two-exciton level. Yet, 
our work extends beyond previous results obtained for single molecules coupled to an MNP, 
because a much more detailed study is carried out than before, considering more different 
conditions. In doing so, we show that the shape of the absorption spectrum in general 
strongly depends on the conditions, in particular, on the detuning between the one-exciton 
resonance and the plasmon resonance and on the ratio between the effective self-interaction 
parameter of the dimer and the dephasing rate of its optical transitions. In the weak-held 
limit, where only the one-exciton transitions play a role, the shape generally is dispersive 
and for a range of exciton-plasmon detunings is Fano-like; close to zero detuning, this shape 
disappears and a deep absorption dip may occur in the broad plasmon absorption line. In 
the strong-held limit, the dimer easily saturates and the Fano-like shape does not occur. 
These results show the sensitivity of the exciton-plasmon hybridization to precise conditions 
and accentuate the possibility to tailor the absorption line shape through appropriate system 
preparation. 

This paper is organized as follows. In the next section, we present the model and the 
theoretical background to treat the optical response of a dimer-MNP hybrid. In Sec. IIIIl 
we discuss the held dependent optical response both of an isolated dimer and a dimer-MNP 
composite (bistability and hysteresis). In Sec. IIVI the Fano-like absorption spectrum of the 
hybrid is discussed, both in the limit of weak and high applied held intensity. We summarise 
in Sec. 0 
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II. MODEL AND THEORETICAL BACKGROUND 


We consider the optical response of a hybrid consisting of a dimer of quantum emitters 
(two level molecules) coupled to a nearby spherical MNP (Ag) embedded in a dielectric 
(dispersionless) background with permittivity Sb- The entire system is subjected to an 
external held of amplitude Eq and frequency uj, £o(t) = Eo(t) cos{ut), oriented along the 
vector d between the centers of the MNP and the dimer. Figure [T](a) schematically shows 
the conhguration. The radius r of the MNP and the MNP-dimer spacing, d, are assumed 
to be small compared to the optical wavelength, thus allowing the use of the quasi-static 
approximation and the point dipole approximation for both particles^*^. Also, the effects of 
retardation are negligible: the interaction between the the MNP and the dimer is dominated 
by the near-held dipole dipole coupling. For the radii of the MNP we consider (r > 10 nm), 
ehects of size quantization of the plasmon spectrum is negligible^ (see Appendix A). Finally, 
we assume that there is no charge exchange between the dimer and the MNP; for d — r < 2 
nm this might occur due to orbital overlap.— 

Like most of the previous work on these types of nanohybrids, we adopt the semiclassical 
approximation, i.e., we treat the optical response of the molecules quantum mechanically, 
while the response of the plasmons in the MNP is treated classically, using its polarizability. 
A full quantum mechanical treatment of the many-electron system that comprises the MNP 
is beyond the scope of this paper. Quantum models treating the MNP essentially with 
a few quantum degrees of freedom recently have been studied22.i2^; while such models are 
interesting, they do not describe effects such as bistable optical response. Thus, we describe 
the MNP by means of its frequency dependent polarizability: 


a{u) = 


^m(^) 

+ 2£f, 


( 1 ) 


where em{^) stands for the dielectric function of the metal. The frequency at which the real 
part of a{u) is minimal determines the plasmon resonance oogp (Frohlich condition). Thus, 
the electronic states of the MNP can be described by a ground state and a continuum of 
excited states with a well dehned plasmon resonance ojgp [See Fig. [11(b), left panel]. For 
the size of our interest corrections to a{u) due to the depolarization shift and radiative 
damping^ are negligible. Likewise, chemical interface damping which important for radii 
below 10 nm can safely be neglected.— 
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FIG. 1. (a) Schematics of a hybrid system of a molecular dimer at a distance d = |d| from a MNP 
subject to an applied field Eq cos uit polarized along d. The dimer consists of two identical quantum 
emitters with equal transition energies Hujq and transition dipole moments and coupling strength 
J. The system is embedded in a homogeneous dielectric host with permittivity £b- (b) Diagram of 
the continuum of energy levels of the MNP (left) and the collective (excitonic) levels of the dimer 
(right). 

We model the dimer as two identical two-level quantum emitters, labeled j = 1,2, with 
ground and excited states indicated as \jg) and \je), respectively; the emitters both have 
a resonance frequency cjq and equally oriented transition dipole moments /x which are, for 
dehniteness, chosen perpendicular to the dimer axis (see Fig. [T]); they have no static dipole 
moment. These two emitters interact resonantly with each other with coupling strength 
J > 0, resulting in four collective (excitonic) eigenstates: the ground state \g) = \lg,2g), an 


5 













antisymmetric and a symmetric one-exciton state, defined through |a) = l/(\/ 2 )(|l 5 ', 2e) — 
\le,2g)) and |s) = 1/ {■\/2){\lg,2e) + \le,2g)), respectively, and the doubly excited state 
|e) = |le,2e). The energies of these states are 0, fvjja = — J, + J) and 

fiLOe = 2huQ, respectively [see Fig. [I][b), right panel]. For the chosen geometry of the dimer, 
indeed the dipople-dipole coupling J is positive, so that the symmetric state is the higher 
state of the one-exciton doublet. The antisymmetric state |a) is dark (zero transition dipole 
moment to the ground state) and does not contribute to the system’s optical response, 
whereas the symmetric one has transition dipole moment D = \/2/x. The transition between 
the doubly excited state |e) and the symmetric one also has the transition dipole moment 
D. Transitions between the dimer’s ground state and doubly excited state are one-photon 
forbiddden and can only be achieved by absorption or emission of two photons. Thus, 
optically, this symmetric dimer represents a three-level ladder-like quantum system. 

We describe the optical dynamics of the dimer within the framework of the density 
matrix formalism. As all dipoles and electric helds are parallel to each other, they can be 
considered as scalars. Then, the set of equations for the density matrix elements pgg, pss, 
Pee, Psg = P*gs, Pes = P*se, ^nd p^g = p*g^ reads (see also Ref. l46|) 

.D£d 


Pgg ‘^'1 Pss + * t i^Psg Psg) , 


h 

pss ‘^'l^es ~\~ i T (Pes Pes T Psg Psg) , 
Pee 2,'^Pee T ^ . \Pes Pes) , 


h 


D£ 


Psg ip^sg T F) Psg “1“ 2“^Pes ~\~ i ^ (Psff Pss T Peg) 


D£ 


pes (P'^es T 2F) pes f ^ i^^es T Peg) , 

■ — (■ _L F'i _L , X 

peg Kp^eg T F) Peg p~ I ^ \Psg Pes) ■ 


( 2 a) 

( 2 b) 

( 2 c) 

( 2 d) 

( 2 e) 

( 2 f) 


Here, 7 is the population radiative rate of a single emitter. Rigorously speaking, 7 should be 
renormalized (enhanced) to (1 -|- F)p, where F is the Purcell enhancement factor due to the 
presence of the closely spaced MNP. In Ref. 1^, it has been shown that F is proportional to 
the real part of the polarizability a{u). Here, the actual frequency range we are interested in 
is around the MNP plasmon resonance, where the real part of a{uj) tends to zero. Therefore, 
the Purcell effect is not important in our case, allowing us to take for the radiative rate its 
free-space value. P = 7 -I-P', with P' denoting the pure dephasing rate of the coherences (off- 
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diagonal density matrix elements); ojgg, ojes, and ojeg are the frequencies of the corresponding 
transitions: \s) -)■ l^f), |e) ^ \s), and |e) -)■ l^f); Zsg = Pss - pgg and = pee - Pss 
are the population differences between the states indicated in the subscripts. Furthermore, 
Sd denotes the electric field acting on the dimer; this field is the sum of the external field 
So = Eocos{ut) and the near held produced by the MNP at the position of the dimer: 


V 


where eo is the vacuum dielectric constant and Vmnp is the MNP’s dipole moment induced 
by the external held Sq and the held produced by the dimer at the position of the MNP: 


'Pmnp — -^^o^bOiiuj) 


Eoe-^^^ + 


D 


r (^Psg E Pes 


+ c.c. . 


(4) 


TieoSbdS’ 

Here the second term in square brackets comes from the held produced by the dimer dipole 
at the position of the MNP: 

Vd = D{psg + Pes + C.C.) (5) 

From Eqs. ([2]) - dS]), we obtain the following expression for the held acting on the dimer: 


Sn — 


1 + 


a{u) 

2n(E 


Eoe 


—iu)t 


+ 


a{u)D 


\Psg “1“ Pes) E C.C. 


( 6 ) 


The term in square brackets in Eq. ([6]) describes the renormalization of the external held 
amplitude Eq due to the presence of the MNP, while the second one is the self-action of the 
dimer via the MNP. As will be shown below, the latter may drastically modify the optical 
response of the composite compared to its separate components. 

For the case of cu ~ cuq and the exciton splitting 2J/h cu, cuq, the rotating wave 
approximation can safely be used to simplify the treatment. To this end, we substitute in 
Eqs. fl2a)) - ([2fl) the oh-diagonal density matrix elements Psg,Pes, Peg) and the held Sd in the 
form 


Psg = Rsge-^^^ , Pes = Rese-^‘^^ , Peg = RegC-^^^^ 

Sd = 2^^^ ’ 


with 


En = 


1 + 


2Tid? 


„ a{uj)D , 


(7a) 

(7b) 


(7c) 
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where the amplitudes Rsg^ Res-, Reg-, and Ed are complex functions slowly varying on the 
time scale of After neglecting all terms oscillating with double frequency 2a;, the 

equations for the slow variables read: 


pgg = 2'ypss + i {^*DRsg - ^DR*sg) , (8a) 

Pss = 2-^Z^s + i{^*DRes-^DR*es + ^DR*sg-^*DRsg) , (8b) 

Pee = -27Pee + * {^dRIs “ ^D^es) , (8c) 

Rsg = ~ (r — i^sg) Rsg + ‘^iRes + * {^^jReg ~ ^D^gg) , (8d) 

Res = - (2r - iAes) Res “ i {^oZes + ^oReg) , (8e) 

Reg = - (r - iAeg) Reg + iVLo^Rgg “ Res) ■ (8f) 


Here, Agg = u — Ugg, Aes = oo — ooes, and Aeg = 2u — u)eg = 2{u) — uo) are the detunings out of 
resonance of the one-photon |s) f-)- \g) and |e) f-)- |s), and two-photon |e) -H- \g) transitions, 
respectively. Furthermore, Qd = DED/{2h) is given by 

Hd = ^0 + G{Rsg + Res) ( 9 ) 


with 


flo — 
G = 


\ , ^i^) ' 

2nd\ 

D^a{uj) 

A'K'^heQShcR 


rio, 
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(10a) 

(10b) 


where Qq = DEQ/{2h) = pE^/{\/2K) is the Rabi frequency of the external field. The 
complex-valued constant G = Gr + iGj is the feed-back parameter arising from the dimer- 
dimer self-action via the MNP. It contains all details of the dimer-MNP interaction, material 
properties, and geometry of constituents. 

It is useful to realize that hG = D^a{u)/{ATi^eoebcR) looks like the inductive dipole-dipole 
interaction between the dimer and the MNP, even if it scales with the dimer-MNP distance d 
as d~^] we use here the term ’’inductive” because an MNP, as a classical object, has neither 
a static nor a transition dipole moment. Up to coefficients of order unity and signs, the 
field produced by the dimer’s transition dipole D at the position of the MNP is D/d^. This 
field induces a dipole Da/d? in the MNP, which in turn has an interaction D'^a/d? with the 
dimer’s transition dipole, confirming our interpretation of G. 






Substituting Eq. (|9]) and Eqs. fllOap - (llObp into Eqs. dHal) - ([Sj), we obtain the equations 
of motion for the density matrix elements of the hybrid in the form 


pgg = 2'ypss + ii^lRsg - ^oKg) 

+ i{G*R*gRsg — GResRgg) + “^GjRsgRgg , 


Pss “t" % 


ai(R.. - R.,) + - R’J 

+ 2iGii{R*gRes — RsgR*es) + ‘^Gi{ResR*es ~ RsgR*sg) ) 

Pee = -2-fPee " 2GlResRl + K^oR^ " ^*oRes) 

+ i{GRsgR*es — G*R*gRes) , 

Rsg = ~ [r — GjZgg + HGjiZsg — A^g)] Rgg + 2'yRes + OVL^Reg — Vt^Zgg) 


+ i [G*{Rlg + Rig)Reg — GZggRes] , 

Res = — [2r — GjZes + OGllZes — Ags)] Res ~ iiS^oZes + ^^Reg) 

— i \GZesRsg + G*{Rlg + Rig)Reg] , 

Reg (h l^eg)Reg R '^^o(^Rsg Res) R '^GiyRggRsg ResRes) •> 


(lla) 

(llb) 

(llc) 

(lid) 

(lie) 
(Ilf) 


As follows from Eqs. flllap - flllfl) . the dimer’s self-action, governed by the feedback 
parameter G, gives rise to additional nonlinearities as compared to an isolated dimer. Two 
of these that should be mentioned in particular, are (i) - renormalization of the dimer’s 
transition frequencies, Ugg UsgRGRZgg and ujes —^ oJesRGuZes, and (ii) - renormalization 
of the damping rates of the off-diagonal density matrix elements, E ^ E — GjZgg and 
2 r —)■ 2r — GjZes, both depending on the corresponding population differences. As will 
be shown below, these two play an important role in understanding the hybrid’s optical 
response. The nonlinearities introduced by these renormalizations are similar to those found 
for a two-level system close to a MNP,->2i'i^ except that more levels are involved here. 


III. RESPONSE AT A FIXED FREQUENCY 

In this section, we investigate the optical response of the nanohybrid at a hxed frequency 
and are interested in the dependence on the magnitude of the driving held (i.e., the Rabi 
frequency Qq)- To this end, we numerically solve Eqs. flllap - flllfp . In our calculations, 
we consider a spherical silver MNP of size r = 11 nm and use the dielectric function 
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for Ag as obtained from experimental data.— For the host’s dielectric constant we assume 
Eb = 1.80. Using these data and applying Eq. o. the peak of the plasmon resonance 
is found around husp = 3.23 eV. Figure H] illustrates the full frequency dependence of the 
polarizability a{u}) for this set parameters. For the two-level systems that make up the dimer, 
we assume a resonance frequency hwQ = 3.13 eV, a transition dipole moment p = 0.5 e • nm 
(with e the magnitude of the electron charge), population and coherence relaxation rates of 
7 = 1.43 ns“^ {h'y = 0.94 x 10“^ meV) and F' = 1.54 ns“^ (hF' = 10“^ meV), respectively, 
and a mutual coupling of J = 0.1 eV. 



fluj [eV] 


FIG. 2. Frequency dependence of the real (Re[a]) and imaginary (Im[a]) parts of the silver MNP 
polarizability a, calculated according to Eq. ([1]) with the set of parameters given in the text. 

We performed calculations of the dimer populations Pgg, pss, and Pee as well as the absolute 
value of the normalized amplitude of the dimer’s dipole moment \Res + Rsg\, sweeping 
adiabatically the magnitude of the applied held up and down. For comparison, we did this 
for both an isolated dimer and a dimer-MNP hybrid. The applied held frequency u was set in 
resonance with the bright dimer transition \g) |s), i.e., fkj = tvujsg = fvjjo + J = 3.23 eV. 

Note that for the chosen parameters Ugg coincides with the plasmon resonance Ugp- The 
results are presented in Fig. |3] (isolated dimer) and Fig. H] (hybrid) and will be discussed 
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below. 


A. Isolated dimer 

Inspecting Fig. El^a), one observes that the populations of an isolated dimer manifest the 
standard behavior when adiabatically sweeping the applied held strength flp- Let us hrst 
discuss the route of increasing held magnitude (indicated in the hgure by solid arrows). As 
expected, in the weak held regime (ffo 7 ), the dimer is in its ground state, pgg ~ 1 (black 
curve), while the other states are unpopulated, Pss ~ Pee ~ 0 (red and blue curves). When flo 
increases, population starts to transfer from the ground state to the symmetric one-exciton 
state: Pgg decreases, accompanied by the growth of pss- This occurs when flo approaches 
7 . The population of the doubly excited state still stays close to zero, because the applied 
held is not resonant with the |s) —)■ |e) transition. Increasing flo further, saturates the 
\g) |s) transition, which is evident in Fig. |3](a) from plateaus in the held dependences at 

pgg = Pss = 1/2. This intermediate regime is characterized by ffo ^ 7 , but still flo ^ T/h. 
Because of the latter condition, the population of the doubly excited state |e) |s) remains 

close to zero. 

Only when Oq ~ T/h, the held is strong enough to overcome the oh-resonance detuning 
of the |s) \e) transition and the latter state is expected to get populated. In fact, at 

this held magnitude, also the transition \g) \e) takes place via simultaneous absorption 

of two photons without the creation of population in the s-state as intermediate. In our 
case, J/h = 10® 7 and, indeed, just at flo ~ 10®7, the population pee begins to grow, 
depleting the populations of the symmetric and ground states. In the limit of very high 
helds, flo ^ T/h, the transitions between all states are equally probable and the system 
becomes fully saturated: the populations of all states are the same, pgg = pee = Pee = 1/3. 
When sweeping down the driving held (indicated in the hgure by dashed arrows), all the 
populations trace back their paths, not showing hysteresis or other deviations (cf. Sec. IIII B|) . 

Figure ini(b) presents the held dependence of the absolute value of the normalized ampli¬ 
tude of the dimer’s dipole moment \Res + Rsg\, obtained while adiabatically sweeping flo 
up and down, in the same way as above. As is observed, \Res + Rsg\ shows two peaks and 
three regions of values close to zero. These observations have clear explanations that can 
be inferred from the held behavior of the populations. Fig. |3](a). For Hq ^ 7 the depletion 
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of the dimer’s ground state is minor and the populations of the excited states |s) and |e) 
are small. Hence, also the coherences Rgg and Res are small, the last one in particular. The 
first maximum of |i?es + Rsg\ is reached at Hq ~ 7) i-®-) at the held magnitude for which 
the population transfer from l^f) —)■ |s) is developed well already. This causes the coherence 
Rgg to grow in magnitude and arrive at its maximum value for Hq ~ 7- At these held 
magnitudes, the coherence Res is still almost zero, because the |s) \e) transition does not 

occur yet. In other words, the low-held peak of \Res + Rsg\ is mainly due to the coherence 
Rsg- 

When approaching the intermediate saturation regime {'j ^ Qq J/h), the populations 
Pss and pgg tend to the same value of 1/2. Because the coherence Rgg is proportional to 
Pss — Pgg, it drops to zero in this regime; as Res still is negligible, also \Res + Rsg\ drops to 
zero. In the strong saturation regime, when Hq > J/K fhe transition |s) ^ |e) occurs and 
the coherence Res starts to grow efhciently, giving rise to the second peak of \Res + Rsg\- It 
should be noticed that the coherence Rgg also is not zero now, which explains the diherence 
in amplitudes of the two peaks. Upon further increasing the driving held magnitude, all 
transitions of the dimer approach the truly saturated regime, whence both Rgg and Res 
decrease, and thus \Res + Rsg\ drops to zero as well. Again, sweeping back, no hysteresis is 
observed. 


B. Dimer-MNP hybrid 


The optical response of a dimer-MNP hybrid is shown in Fig. IH Here, we set the center- 
to-center distance between the MNP and the dimer at d = 18 nm; the other parameters are 
as in Sec. IIII A1 Comparing this hgure to Fig. [3l it is obvious that the hybridization leads to 
strong changes of the response. In contrast to the isolated dimer, the response of the hybrid 
depends on the history of sweeping the applied held (up or down): in a certain interval of 
Hq) the quantities Pgg, pss, and \Res + Rsg\ exhibit a hysteresis loop, related to the \g) |s) 

transition. This indicates bistability of the hybrid.— The effect results from the self-action 
of the dimer via the MNP, described by the feed-back parameter G [Eq. (llObll ] and also 
occurs for a hybrid of a MNP and a monomer-*!^, as it only requires one exciton transition. 

The real and imaginary parts of G (denoted G^ and Rj, respectively), give rise to different 


mechanisms for bistability, as has been analyzed in detail in Ref. 


121. In order to achieve 
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bistable response, the values of Gr and Gj should meet certain threshold conditions.—*^ For 
the parameters chosen in Fig. 01 we have G = (27.47 + 936.64 z)F [hG = (0.05 + 1.83 z)meV], 
and indeed the criteria for bistability are met. 

From comparing Figs. ([3]) (a) and (l4])(a), one observes that the field magnitudes at which 
the two saturation steps occur differ. For the hybrid the saturation of the one-exciton 
transition, indicated by the plateau at pee = Pss = 1/2 occurs at a value of flo that is 
larger than for the isolated dimer. The explanation lies in the (population dependent) 
renormalization of the resonance frequency and the coherence relaxation rate due to the 
dimer-MNP coupling: ojgg Ugg + GRZ^g and F ^ F — GjZgg. At low values of hlo, fhe 
population difference Zgg = —1, so that the renormalization brings the effective resonance 
of the dimer away from the frequency of the applied held and at the same time increases 
the effect of damping. Realizing that both Gr and Gj are much larger than F, these 
renormalizations imply that a higher value of Hq is required to saturate the transition. 
Oppositely, for the hybrid, the two-exciton transition saturates at a lower intensity than for 
the isolated dimer. The explanation lies in the enhancement of the held acting on the dimer 
due to the exciton plasmon interaction [see Eq. ([9])]. All the above ehects also manifest 
themselves in shifts of the peaks in the held dependencies of \Res + Rsg\ [Figs. ([3])(b) and 

0 (b)]. 

To get more insight into the held dependence of the hybrid’s response, we calculated the 
average number of excited emitters, given by n = pss + ‘2pee, for three diherent frequencies of 
the applied held: hu = hugg = huo + J = 3.23 eV (in resonance with the \g) —)> |s) transition 
as well as the MNP plasmon Ugp = 3.23 eV), hw = hojQ = 3.13 eV (in resonance with the 
two-photon transition l^f) |e)), and hw = hojes = hcuo — J = 3.03 eV (in resonance with 
the |s) \e) transition). Figure 0] shows the results for an isolated dimer [G = 0) and a 

dimer-MNP hybrid (G 7 ^ 0), adiabatically sweeping the applied held magnitude hlo up and 
down. 

For huj = tvjJsg = 3.23 eV [Fig[5][a)], n clearly displays bistable behavior in the hybrid. 
Upon sweeping Hq up, initially n remains close to zero, meaning that the dimer is in its 
ground state \g). At some critical value of ffo ~ the \g) |s) transition, dissimilar 

to an isolated dimer, is abruptly saturated (n = 1 / 2 ), maintaining this state within some 
interval of flo- Upon increasing flo beyond that interval, all dimer transitions become sat¬ 
urated. Sweeping Hq down again, the transition from the fully saturated regime (n = 1) 
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to the intermediate one (n = 1/2) shows no hysteresis, while the second step down to the 
unsaturated regime {n = 0) again occurs abruptly, but at a lower value of flo than the 
upward step. Thus, the average number of excited emitters exhibits hysteresis, conhrming 
that the dimer coupled to the MNP is a bistable system, where the bistability occurs only 
in the \g) |s) channel and does not appear in the |s) \e) channel. No multistable^ 

behaviour is found. 

If the driving held is tuned away from the \g) —?■ |s) resonance. Fig. [5])b), no bistability 
occurs in the hybrid’s optical response. However, the effect of local-held enhancement due 
to the plasmon excitation is seen for all cases: the applied held magnitude Hq required to 
completely saturate the system always reduces signihcantly in the presence of the MNP. We 
hnally note that for the held frequency ho; = hwo = 3.13eV, the system gets fully saturated at 
the lowest intensity as compared to all others excitation conditions. This can be understood 
from the fact that then the applied held is in exact resonance with the \g) \e) two-photon 

transition. 


IV. HYBRID ABSORPTION SPECTRA 

In this section, we analyse the absorption spectrum of the dimer-MNP hybrid and show 
that its shape strongly depends on the detuning between the one-exciton resonance ujgg and 
the plasmon resonance as well as the ratio |G|/P. We will distinguish between the 
weak-held limit Oq ~ 7, where the two-exciton state |e) is irrelevant, and strong-held limit 
(Oq ^ 10 ^ 7 ), where the state |e) participates in the process. 

The spectra are computed from the time average of the power (energy per second) Q 
absorbed by the hybrid; Q consists of two contributions, one from the dimer (Qdim) and 
the other from the MNP (Qmnp): Q = Qdim + Qmnp- From now on, we will refer to the 
frequency dependence of all Q’s as the absorption spectra of the corresponding parts. 

The absorption of the dimer results from the creation of excitons, while in the case of the 
MNP, it is due to charge oscillations (plasmons). Here, we do not include the contribution 
of scattering by the MNP, because for small particles (r << A) this ehect is negligible. The 
power absorbed by the dimer (MNP) is given by:—: 

QdIM(MNP) = 2^0l^fh)IM(MNP)Ini[XDIM(MNP)] IF'cT ) (12) 
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where Xdim(mnp) is the susceptibility of the dimer (MNP), defined as the proportionality 
coefficient between the polarization (dipole moment per unit volume) and the applied field 
amplitude Eq, and Pdim(mnp) is the volume of the dimer (MNP). 

The dimer susceptibility Xdim is calculated quantum mechanically and has the form 


_ D{Res + Rsg) 

^ 2Vl>iMeo-Bo 


( 13 ) 


Note that the volume Vdim attributed to the dimer is somewhat arbitrary. However, in the 
absorbed power, Eq. flT^ . this volume drops out so that this arbitrariness is irrelevant. 

The MNP susceptibility Xmnp is computed within the framework of classical electrody¬ 


namics: 


Xmnp — 


EfjCy 


1 + 


D 


-{Re 


R 


sg) 


(14) 


VmNP L 

The second term in square brackets represents a correction to the external held Eq by the 
held produced by the dimer at the position of the MNP. Having obtained Xdim and Xmnp, 
we can now calculate the absorption spectrum Q of the dimer-MNP composite. 


A. Weak field regime 

First, we address the linear absorption spectrum of a dimer-MNP hybrid under weak held 
excitation. In particular, we set the Rabi frequency Hq = 7 - Then, the |s) \e) transition 

is not involved in the hybrid’s absorption. We consider diherent positions of the frequency of 
the dimer’s one-exciton resonance ujsg with respect to the plasmon resonance frequency ujsp 
and vary the dephasing rate P'. The other parameters used in the calculations are identical 
to those described above. 

Figure [ 6 ] shows the absorption spectra obtained for a dimer-MNP hybrid with center- 
to-center distance d = 18 nm (strong coupling regime) at three values of the one-exciton 
resonance: hwag = 3.14 eV (panel (a), below the plasmon resonance), fvjJsg = 3.23 eV (panel 
(c), at the plasmon resonance), and hjJsg = 3.32 eV (panel (e), above the plasmon resonance). 
In each case, results are shown calculated for three different dephasing rates: P' = 1.08 7 
{hV = 10“^ meV), P' = 1.08x 10^ 7 {hV = 0.1 meV), and P' = 1.08x 10^ 7 {HT' = 1 meV). 
In the absence of the exciton-plasmon interaction, the spectrum would be a simple Lorentzian 
lineshape centered at = 0 and with a width given by 7 -f P', superimposed on the broad 
Lorentzian plasmon peak centered at Agg = Usp—oOsg- As is clearly seen, the exciton-plasmon 
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interaction, given by the coupling parameter G, strongly affects the spectrum, which now 


exhibits an asymmetric dispersive shape with a dip to almost zero absorption (for small F' 
values), resembling the well-known Fano resonance^. Recall that the essence of the original 


quantum Fano effect - response of a system comprising a discrete level interacting with a 
dark quantum continuum - consists of the quantum interference of the mixing amplitudes 
from the bare states to the hybridized (due to the interaction) collective states. In our case, 
the continuum is a classical system with losses: surface plasmon of the MNP. Here, the 


quantum picture can not be directly applied, but by analogy, the interference of the applied 
field and the field produced by the hybrid leads to the Fano-like spectrum.->^ Whether the 


interference is constructive (peak) or destructive (dip) depends on the detuning away from 
the exciton resonance. 

Panels (b), (d), and (f) in Figj 6 ]show the frequency dependence of the coupling parameter 
G in the relevant frequency ranges. Evidently, when the exciton resonance lies in one of 
the tails of the plasmon peak, both G^ and Gj are almost constant [Figs. | 6 ][b) and (f)]. 
However, when the exciton is (close to) resonant with the plasmon, G (in particular Gr) 
shows a pronounced frequency dependence [Fig. | 6 ](d)]. Note that, when Usg lies in the low- 
frequency tail of the plasmon peak, the real part of a{u) is positive and therefore Gr is 
positive as well [Figj 6 ][b)]. The peak in the spectrum then occurs at lower frequency then 
the dip [Figj 6 ][a)]. For cjsg in the higher frequency tail of the plasmon peak, Gr is negative 
and the peak in the spectrum occurs at higher frequency than the dip [Figs. [ 6 ][e) and (f)]. 

We will Erst make several observations regarding the coupled exciton-plasmon system 
under resonant conditions. Whether or not the dimer-MNP system is coherently coupled, 
depends, in fact, on the ratio of the dimer-MNP coupling strength |G| and the dephasing rate 
P = 7 -I-P'. To investigate this, consider the situation where the dimer resonance approaches 
the maximum of the plasmon peak {usg —t (u^p), for which the spectrum shows only a dip in 
the absorption and no peak [Fig. EJ^c)]. In this case, the imaginary part Im[a] acquires its 
maximum, whereas Re [a] tends to zero [see Fig. [2]. The latter allows one to derive a simple 
analytical formula for the susceptibility of the dimer and the MNP (see Appendix |B]) which 
leads to the totally absorbed power by the hybrid in the form 



( 16 ) 


Here, ai = Im[Q;]. The first and second terms in brackets correspond to the absorption of 
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the MNP and the dimer, respectively. As is seen from Eq. ffTSj) . the minimum value of the 
hybrid absorption (reached in the dip) is governed by the ratio Gj/T. 

The physical explanation is as follows. In the strong coupling regime [Gj S> P), both 
the dimer and the MNP acquire a large collective damping, given by Gi, thus reducing 
the hybrid’s absorption. The dephasing counteracts the collective nature of the hybrid’s 
states and works towards recovering the absorption near which finally, at Gi P, is 
represented by the absorption of the isolated MNP and dimer (first and second terms in the 
parentheses of Eq. [151 respectively). This is the origin of the dependence of the dip depth 
on the ratio G//P, reflected in Fig. | 6 ](c). 

The role of the ratio |G|/P may also be investigated by varying G. Figure [7] illustrates 
the crossover from the weak to the strong coupling regime upon changing the dimer-MNP 
distance d from large to small values. In the calculations, we used the same set of parameters 
as in Fig[ 6 ](e) (we are thus in the high-energy plasmon tail), except that the pure dephasing 
rate was hxed at P' = 1.08 7 {hT' = 10“^ meV) and the dimer-MNP distance was varied 
from d = 38 to 18 nm. As is seen from Fig. [71 at large dimer-MNP distances, the absorption 
spectrum of the dimer-MNP composite displays a narrow peak (bright red) which coincides 
with the resonance of an isolated dimer (shown by the vertical white dashed line), and its 
width is just P'. 

Upon decreasing d, the dimer peak weakens, because the dimer excitation starts to in¬ 
teract with the MNP and to share its energy with the plasmons. Around approximately 
d = 25 nm, the crossover to the strong coupling regime occurs and the system falls into the 
Fano regime. This is reflected in the absorption spectrum: a peak (red dotted line) and a 
dip (white dotted line) are clearly visible. Note also that upon decreasing d, the peak is 
blue-shifted with respect to the isolated dimer resonance and broadens. The black dashed 
curve highlights the d-dependence of the frequency shift. 

Within a semiclassical picture (a quantum emitter coupled to a classical MNP), the shift 
and broadening have a natural explanation. Recall that the dimer resonance frequency Ugg 
and the dephasing rate P are renormalized due to the exciton-plasmon coupling: Ugg —)■ 
Usg + GnZsg and P ^ P — GjZgg. Under the conditions considered here (weak excitation 
and blue-tail location of the dimer resonance with respect to the plasmon peak), both Gr 
and Z = —1 are negative, while Gj is always positive [see plots (b), (d), and (f) in Fig. [ 6 ]. 
This indeed predicts a blue-shifted renormalized frequency and an increase of the line width; 
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moreover, both effects increase when d decreases. 


B. Strong field regime 

In this section, we address the spectrum of the dimer-MNP hybrid in the limit of a strong 
electric field, when the two-photon transition from l^f) \e) comes into play. The hybrid’s 
parameters were chosen as in Sec. IIV Al with a fixed center-to-center distance d = 18 nm 
and restricting to the dimer’s dephasing rate T' = 1.08 7 {hT' = 10“^ meV). 

Figure [H] displays the absorption spectrum of the hybrid calculated for two values of the 
external field magnitude: hlo = (blue solid curve) and hlo = lO^'^y (black dashed 

curve). For hlo = lO^y, the system still exhibits a dip in the absorption when the frequency 
Ugg of the isolated dimer’s \g) |s) transition coincides with the MNP plasmon resonance 

Usp [Fig. [H] (a)]. However, the dip is much less pronounced than in the weak field limit (see 
Fig. E]), and does not drop to zero absorption anymore. This occurs due to saturation of the 
dimer’s \g) —?> |s) transition, which also reduces the dimer-MNP interaction. The dip is also 
visible for cjsg > oogp [see Fig. [ 8 ] (b)]. 

Upon increasing the applied field magnitude to flo = lO^'^y, the dip in the hybrid’s 
absorption is completely smeared, independently of the excitation frequency. This results 
from an even stronger saturation of the dimer’s \g) |s) transition, which also experiences 

power broadening. At these strong fields, an additional peak in the isolated dimer spectra 
appears at the frequency of the two-photon transition JtjjJq, which equals, respectively, 3.13 
eV [FigJH] (a), right insert] and 3.22 eV [Figj 8 ][b), right insert]. This additional peak does 
not produce any discernible feature in the hybrid’s absorption spectrum, meaning that the 
direct two-photon transition does not interact strongly with the MNP plasmon resonance. 
This results from the fact that this transition is only allowed in second-order perturbation 
theory, leading to a small effective transition dipole. 


V. SUMMARY 

In this work, we theoretically studied the optical response of a nanohybrid comprising 
a symmetric quantum dimer in close proximity to a MNP. The symmetric dimer consists 
of two identical quantum emitters (two-level molecules) and acts as an optical three-level 
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system, with a ground state \g) (both emitters in the ground state), an optically active one- 
exciton state |s) (the symmetric linear combination of states in which one emitter is excited 
and the other one is in the ground state), and a doubly excited state |e) (both emitters 
excited). The antisymmetric one-exciton state is optically dark and may be ignored. The 
doubly excited state is optically active, however, it only can be excited via a two photon 
process, either through the one-exciton state as intermediate or by the direct absorption of 
two photons. We derived equations for the dimer’s 3x3 density matrix describing its optical 
dynamics in the presence of a metal nanoparticle and analyzed these equations to study the 
optical properties of the hybrid. 

We found that in a certain range of system parameters, the hybrid may exhibit bistabil¬ 
ity, but that this effect only occurs for the one-exciton transition, not for the two-exciton 
transition. The mechanism of this effect is identical to that of a single quantum emit¬ 
ter close to an MNP.— The two-exciton transition exhibits excitation enhancement which 
causes saturation of the hybrid system at held intensities signihcantly lower than for an 
isolated dimer. The later indicates enhancement of the two photon absorption as has been 
observed experimentally.— No multi-stability occurs. 

Another fascinating property of the dimer-MNP hybrid is its Fano-like absorption spec¬ 
trum, with a shape that strongly depends on the spectral location of the dimer’s one-exciton 
absorption peak with respect to the MNP’s plasmon resonance and the ratio of coupling pa¬ 
rameter and the dephasing rate. In particular, its basic shape depends on whether the 
exciton resonance lies at the low or the high frequency side of the plasmon resonance, or 
just at its maximum. In all cases the spectrum shows almost zero absorption at some fre¬ 
quency, making the hybrid transparent at this point. This effect is especially pronounced at 
the maximum of the MNP plasmon resonance. Here, the hybrid absorption shows a narrow 
dip to almost exactly zero. 

Upon increasing the applied field intensity, the dip becomes less pronounced and disap¬ 
pears when the hybrid gets saturated. Similar effects occur if the one-exciton energy lies in 
one of the wings of the plasmon peak. In addition, the shape of the lineshape is sensitive to 
the coupling between the dimer’s exciton transition and the plasmons. All this allows one 
to tailor the Fano-like shape of the hybrid’s absorption spectrum by controlling the spectral 
and spatial separation of both constituents and by changing the input power. 

We hnally note that in this paper, we considered the limit where the magnitude of the 
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coupling parameter G is smaller than the spectral detuning between the \g) |s) and 

the |s) \e) transitions of the quantum dimer. Preliminary calculations show that in the 

opposite limit, the hybrid shows a rich nonlinear dynamics, including auto-oscillations and 
chaos. This will be the subject of a forthcoming article. 
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Appendix A: Size quantization of plasmons in Ag MNPs 

Here, we address the question of size quantization of plasmons confined in a small volume. 
The importance of size effects for an MNP can be estimated from comparing the energy level 
spacing in the vicinity of the Fermi energy Sp with the level width Pp and thermal energy 
kpT {kp is the Boltzman constant). The spacing is given by^ 

Ae = {2/3)eF/N , (Al) 

where N is the number of electrons in a nanoparticle. For an Ag nanoparticle, Sp = 5.48 eV 
and the electron density N/Vmnp = 5.85 x 10^® m“^, where Vmnp is the Ag nanoparticle 
volume. Then, for a spherical Ag nanoparticle of size r = 11 nm (our case), we get N = 
3.26 X 10® and, respectively, for the energy spacing Ae 0.01 meV. The thermal energy 
kpT ^ 25 meV (at room temperature), whereas Fp 0.1 fs ^ 65 meV.— Thus, both 

kpT and Fp are much larger than Ae, making the size quantization of conhned plasmons 
irrelevant in our case. In a recently published paper,— it has been shown experimentally 
that size effects in Ag nanoparticles only play a role for radii smaller than to approximately 
5 nm. 
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Appendix B: Susceptibilities xdim (dimer) and xmnp (MNP) at uigg ojgp in the 
weak-field regime 


In the weak field regime (fl^i ^ 721 ); the |s) \e) transition is not involved in the dimer- 

MNP hybrid absorption. Moreover, the ground state depletion in this limit is negligible, 
meaning that pgg ~ 1, pss ~ 0. Thus, the only relevant equation in this case is Eq. fllldj) for 
the coherence Rsg, which now takes the form 

Rgg = — [r + G/ — i{^sg + Rsg — *h2o • (Bl) 


We are interested in the steady-state solution to this equation, which reads 


^ r -|- Gi — i(Asg -|- Gpi) 


1 + 


a 


2TTd? 


a 


0 , 


(B2) 


where we used fio = [1 + a/{2 'K(R)]VLo. Having obtained the analytical solution for Rgg, one 
can write down the expressions for the dimer and MNP susceptibilities xt> and Xmnp, the 
imaginary parts of which determine the absorption spectra. 

First, we calculate Xmnp- From Eqs. flTTl) and flB2ll . it is obtained as 


Xmnp — 
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MNP 


1 + 
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1 + 
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271(P 


(B3) 


The general expression for Im[xMNp] is bulky and hard to analyse analytically. However, at 
the maximum of the MNP plasmon peak, u = Ugp, and assuming perfect resonance between 
the dimer and the MNP plasmon, Ugg = Ugp, one can derive a simple formula for Im[Q!]. 
Realizing that under these conditions both Re[a] and Gr are equal to zero [see Fig. [Sl[d)], 
we obtain 


Ini[x 
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where = Im[Q!]. Because the factor that multiplies 1/(P -|- Gj) in the second term within 
square brackets is nothing else than Gj, one finally gets 
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Within the same approximations as above, one obtains for the imaginary part of the 
dimer susceptibility 
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FIG. 3. Field dependent response of an isolated dimer: (a) the populations pgg, pss, and p^e', (b) the 
absolute value of the normalized amplitude of the dimer’s dipole moment \Res + Rsg\- Solid and 
dashed arrows show the direction of adiabatically sweeping the applied field magnitude Oq up and 
down, respectively. The frequency of the field is taken to be in resonance with the dimer’s bright 
1^) —)■ |s) transition, Hlo = Hujsg = hujQ + J = 3.23 eV. 


25 













1 


-► 

' 


' 

(a) 

* — 

A 



- Pgg 

1 


yf 

Pee 





-> 

- 

- 

tk - 

<: - 


-► 


, 


, 

■tt -;■ 





logio(r2o/7) 



logio(no/7) 


FIG. 4. Same as in Fig. [3l but now for a dimer-MNP hybrid with the centre-to-centre distance 
between the MNP and the dimer of d = 18 nm. (a) The populations pgg, pss, and pee- (b) Absolute 
value of the normalized amplitude of dimer’s dipole moment \Res + Rsg\- 
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FIG. 5. Field dependence of the average number of excited emitters, n = pss + ‘^Pee^ for an isolated 
dimer (G = 0) and for a dimer coupled to a MNP (G 7 ^ 0), calculated for three types of excitation: 
(a) in resonance with the \g) —)• |s) transition, Huj = hwgg = hujQ + J = 3.23 eV, (b) in resonance 
with the two-photon transition l^f) —>■ |e), with frequency hw = hwQ = 3.13 eV, and in resonance 
with the |s) —)• |e) transition, hw = — J = 3.03eV. As before, the pure dephasing constant 

equals F' = I.O 87 . Solid and dashed arrows show the direction of adiabatically sweeping the applied 
field magnitude Qq. 27 
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FIG. 6 . Absorption spectra of the dimer-MNP hybrid calculated for different locations of the 
isolated dimer resonance cogg with respect to the MNP plasmon resonance cogp = 3.23 eV. (a) hcogg = 
3.14 eV (below the plasmon resonance), (c) hugg = 3.23 eV (at the plasmon resonance), and 

(e) hujQ = 3.32 eV (above the plasmon resonance). The inserts illustrate these three positions: 

the dashed curves represent the MNP plasmon peak, while the vertical lines indicate the isolated 

dimer resonance. Each panel displays the spectra of the hybrid computed for different dephasing 

rates: P' = 1.08 7 {hT' = 10“^ meV) (black solid curve), P' = 1.08 x 10^ 7 {hT' = 0.1 meV) (red 

dashed curve), and P' = 1.08 x 10^7 {hT' = 1 meV) (blue dash-dotted curve). Panels (b), (d), and 
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(f) show the spectra of the real and imaginary parts of the coupling parameter G. 






















































FIG. 7. Contour plot of the dimer-MNP hybrid absorption spectra Q in the vicinity of the isolated 
dimer resonance Ugg calculated for different center-to-center distances d. Decreasing d corresponds 
to increasing the dimer-MNP coupling. The dashed black curve shows the d-dependence of the 
resonance frequency blue shift Ugg + Gr {Gr > 0). The red and the white dotted lines highlight 
the position of the peak and the dip, respectively. The white dashed vertical line denotes the 
resonance frequency of the isolated dimer uigg- The set of parameters are: hujQ = 3.32 eV (above 
the plasmon resonance), P' = 1.08 7 {hV = 10“^ meV). 
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FIG. 8. Absorption spectra of the dimer-MNP hybrid in the strong field regime calculated for 
two values of the applied field strength: Qq = 10^7 (blue solid curve) and Oq = 10^'^7 (black 
dashed curve). Two dimer resonance frequencies are considered: (a) hugg = 3.23 eV (at the 
plasmon resonance) and (b) fujsg = 3.32 eV (above the plasmon resonance). Inserts show spectra 
of the isolated dimer (solid curve) and of the MNP (dashed curve) for Qq = 10^7 (left insert) and 
Ho = 10^'^7 (right insert). 
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